Phonon-mediated decay of an atom in a surface-induced potential 
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We study phonon-mediated transitions between translational levels of an atom in a surface-induced 
potential. We present a general master equation governing the dynamics of the translational states 
of the atom. In the framework of the Debye model, we derive compact expressions for the rates 
for both upward and downward transitions. Numerical calculations for the transition rates are 
performed for a deep silica-induced potential allowing for a large number of bound levels as well 
as free states of a cesium atom. The total absorption rate is shown to be determined mainly by 
the bound-to-bound transitions for deep bound levels and by bound-to-free transitions for shallow 
bound levels. Moreover, the phonon emission and absorption processes can be orders of magnitude 
larger for deep bound levels as compared to the shallow bound ones. We also study various types of 
transitions from free states. We show that, for thermal atomic cesium with temperature in the range 
from 100 [iK to 400 /iK in the vicinity of a silica surface with temperature of 300 K, the adsorption 
(free-to-bound decay) rate is about two times larger than the heating (free-to-free upward decay) 
rate, while the cooling (free-to-free downward decay) rate is negligible. 

PACS numbers: 34.50.Dy,33.70.Ca 
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I. INTRODUCTION 



Over the past few years, tight confinement of cold 
atoms has drawn considerable attention. The interest in 
this area is motivated not only by the fundamental na- 
ture of the problem, but also by its potential applications 
in atom optics and quantum information. A method 
for microscopic trapping and guiding of individual atoms 
along a nanofiber has been proposed Surface-atom 
quantum electrodynamic effects have constituted another 
interesting area, where a great deal of work has been 
carried out. Modification of spontaneous emission of 
an atom and radiative exchange between two distant 
atoms [H mediated by a nanofiber have been investigated. 
Surface-induced deep potentials have played a major role 
and have received due attention in recent years. Oria 
et al. have studied various theoretical schemes to load 
atoms into such potentials 0, [||. A rigorous theory of 
spontaneous decay of an atom in a surface-induced po- 
tential invoking the density-matrix formalism has been 
developed @ . The role of interference between the emit- 
ted and reflected fields and also the role of transmission 
into the evanescent modes were identified. Further cal- 
culations on the excitation spectrum have been carried 
out Q. Bound-to-bound transitions were shown to lead 
to significant effects like a large red tail of the excita- 
tion spectrum as compared to the weak consequences of 
free-to-bound transitions. A crucial step in this direction 
was the experimental observation of the excitation spec- 
trum and the channeling of the fluorescent photons along 
the nanofiber Q, opening up avenues for novel quantum 
information devices. 

In most of the problems involving surface-atom inter- 
action, the macroscopic surface is usually kept at room 
temperature. Thus the pertinent question that can be 
asked is what would be the effect of heating on the cold 



atoms. It is understood that transfer of heat to the 
trapped atoms will lead to a change in the occupation 
probability of the vibrational levels as well as their co- 
herence. Phonon-induced changes in the populations of 
the vibrational levels have been studied by several groups 
[E H El • In a nice and compact treatment based on the 
dyadic Green function and the Fermi golden rule, Henkel 
et al. showed that the effects can be very different de- 
pending on the nature of the atomic/molecular species 
|9j. The time scales for various species were estimated. 
It should be stressed that the trap considered by Henkel 
et al. was not necessarily a surface trap and misses out 
on many of the aspects of the surface-atom interaction 
[9j. Based on the assumption that the surface-atom in- 
teraction can be represented by a Morse potential, the 
phonon-mediated decay was estimated by Oria et al. [H . 
Their estimate was based on the formalism developed by 
Gortel et al. [10]. However, all the previous theories 
focus on only the transition rates and thus are not gen- 
eral enough. In this paper, we present a general density- 
matrix formalism to calculate the phonon-mediated de- 
cay of populations as well as the changes in coherence. 
We derive the relevant master equation for the density 
matrix of the atom. We emphasize that our density- 
matrix equation describes the full dynamics of the cou- 
pling between trapped atoms and phonons and does not 
assume any particular form of the trapping potential. 
Under the Debye approximation, we derive compact ex- 
pressions for the phonon-mediated decay rates. Numer- 
ical calculations are carried out assuming the potential 
model considered in Q . In contrast to the previous work, 
we include a large number of vibrational levels due to the 
deep surface-atom potential. We show that there can be 
significant differences in the decay rates when the initial 
level is chosen as one of the shallow or deep bound levels. 
We also calculate and analyze the decay rates for various 
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types of transitions from free states. 

The paper is organized as follows. In Sec. |TT] we de- 
scribe the model. In Sec. Mil we derive the basic dynam- 
ical equations for the phonon-mediated decay processes. 
In Sec. IIVI we present the results of numerical calcula- 
tions. Our conclusions are given in Sec. [Vj 



II. DESCRIPTION OF THE MODEL SYSTEM 

We assume the whole space to be divided into two re- 
gions, namely, the half-space x < 0, occupied by a nondis- 
persive nonabsorbing dielectric medium (medium 1), and 
the half-space x > 0, occupied by vacuum (medium 2). 
We examine a single atom moving in the empty half- 
space x > 0. We assume that the atom is in a fixed 
internal state \i) with energy Twji. Without loss of gen- 
erality, we assume that the energy of the internal state 
\i) is zero, i.e. Ui — 0. We describe the interaction be- 
tween the atom and the surface. We first consider the 
surface-induced interaction potential and then add the 
atom-phonon interaction. 



A. Surface-induced interaction potential 



Here m is the mass of the atom. In the numerical ex- 
ample with atomic cesium, we have m = 132.9 a.u. 
= 2.21 x 10 -25 kg. The eigenvalues E v are the center- 
of-mass energies of the translational levels of the atom. 
These eigenvalues are the shifts of the energies of the 
translational levels from the energy of the internal state 
\i). Without loss of generality, we assume that the 
center-of-mass eigenfunctions (f v {x) are real functions, 
i.e. ipl(x) = <p v {x). 

In Fig. [TJ we show the potential U (x) and the wave 
functions tp v {x) of a number of bound levels with en- 
ergies in the range from —1 GHz to —5 MHz. We also 
plot the wave function of a free state with energy of about 
4.25 MHz. In order to have some estimate about the spa- 
tial extent of a wave function ip v {x), we define a crossing 
point cCcross, which corresponds to the rightmost solution 
of the equation U(x) = £ u . Note that, for shallow lev- 
els, the wave function generally peaks close to the point 
Across- We plot the eigenvalue modulus \£ v \ and the cross- 
ing point x cr oss in Figs.[DJa) andHfb), respectively. It is 
clear from the figure that, for v in the range from to 
300, the eigenvalue varies dramatically from about 158 
THz to about 322 kHz, while the wave function extends 
only up to 170 nm. 



In this subsection, we describe the interaction between 
the atom and the surface in the case where thermal vi- 
brations of the surface are absent. The potential en- 
ergy of the surface-atom interaction is a combination of 
a long-range van der Waals attraction and a short-range 
repulsion [11| • Despite a large volume of research on the 
surface-atom interaction, due to the complexity of sur- 
face physics and the lack of data, the actual form of the 
potential is yet to be ascertained [III- For the purpose 
of numerical demonstration of our formalism, we choose 
the following model for the potential 0, [ll| : 



U(x) 



9i 

X 3 



(1) 



Here, C3 is the van der Waals coefficient, while A and 
a determine the height and range, respectively, of the 
surface repulsion. The potential parameters C3, A, and 
a depend on the nature of the dielectric and the atom. 
In numerical calculations, we use the parameters of fused 
silica, for the dielectric, and the parameters of ground- 
state atomic cesium, for the atom. The parameters for 
the interaction between silica and ground-state atomic 
cesium are theoretically estimated to be C3 = 1.56 kHz 
^m 3 , A = 1.6 X 10 18 Hz, and a = 53 mf 1 @. 

We introduce the notation ip u (x) for the eigenfunc- 
tions of the center-of-mass motion of the atom in the 
potential U{x). They are determined by the stationary 
Schrodinger equation 



2m dx 2 



ip v (x) = E v ip v {x). 



(2) 



FIG. 1: Energies and wave functions of the center-of-mass 
motion of an atom in a surface-induced potential. The pa- 
rameters of the potential are C3 = 1.56 kHz fim 3 , A = 
1.6 x 10 18 Hz, and a = 53 nm" 1 . The mass of the atom 
is m = 2.21 x 10 -25 kg. We plot bound levels with energies 
in the range from —1 GHz to —5 MHz and also a free state 
with energy of about 4.25 MHz. 



FIG. 2: Eigenvalue modulus \£ v \ (a) and crossing point x CIOSS 
(b) as functions of the vibrational quantum number v. The 
parameters used are as in Fig. ff] 

We introduce the notation \v) = \<p v ) and uj u — E v j% 
for the state vectors and frequencies of translational lev- 
els. Then, the Hamiltonian of the atom in the surface- 
induced potential can be represented in the diagonal form 



— TlLJ u <7 ulJ . 



(3) 



Here, a vv = |^)(^| is the population operator for the 
translational level v. We emphasize that the summation 
over v includes both the discrete {£ u < 0) and continuous 
(£„ > 0) spectra. The levels v with £„ < are called 
the bound (or vibrational) levels. In such a state, the 
atom is bound to the surface. It is vibrating, or more 
exactly, moving back and forth between the walls formed 
by the van der Waals part and the repulsive part of the 
potential. The levels v with £„ > are called the free (or 
continuum) levels. The center-of-mass wave functions of 
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the bound states are normalized to unity. The center-of- 
mass wave functions of the free states are normalized to 
the delta function of energy. 



B. Atom— phonon interaction 

In this subsection, we incorporate the thermal vibra- 
tions of the solid into the model. Due to the thermal 
effects, the surface of the dielectric vibrates. The surface- 
induced potential for the atom is then U(x — x s ), where 
x s is the displacement of the surface from the mean po- 
sition (x s ) — 0. We approximate the vibrating potential 
U(x — x s ) by expanding it to the first order in x s , 



with 



U(x-x s ) = U(x) - U'(x)x s . 



(4) 



The first term, U (x) , when combined with the kinetic 
energy p 2 /2m, yields the Hamiltonian Ha [see Eq. 
which leads to the formation of translational levels of 
the atom. The second term, —U'(x)x s , accounts for the 
thermal effects in the interaction of the atom with the 
solid. Note that the quantity F = —U'{x) is the force 
of the surface upon the atom. Hence, the force of the 
atom upon the surface is — F = U'{x) and, consequently, 
U'(x)x s is the work required to displace the surface for 
a small distance x s . 

It is well known that, for a smooth surface, the gas 
atom interacts only with the phonons polarized along 
the x direction [10]. In the harmonic approximation, we 
have 



2MNu> c 



1/2 



(6 q e iqR 



6 q e" 



iqR\ 



(5) 



Here, M is the mass of a particle of the solid, N is 
the particle number density, w q and q are the frequency 
and wave vector of the cc-polarized acoustic phonons, re- 
spectively, R = (0, y, z) is the lateral component of the 
position vector (x, y, z) of the atom, and 6 q and 6 q are 
the annihilation and creation phonon operators, respec- 
tively. Without loss of generality, we choose R = 0. 
Meanwhile, the operator U 1 can be decomposed as V = 
12 VV ' a w {v\U'\v'), where a vv > = \v){v'\ is the operator 
for the translational transition v <-> v' . Hence, the en- 
ergy term —U'(x)x s leads to the atom-phonon interac- 
tion Hamiltonian [irj| 



/OJn 



-S(b q 



with 



S — ^ ' Qvy'Vyv' ■ 



(6) 



(7) 



Here we have introduced the atom-phonon coupling co- 
efficients 



(8) 



F vv i= — I ip v (x)U'{x)<p v i{x)dx 



(9) 



being the matrix elements for the force of the surface 
upon the atom. We note that F vv i = —mu>' 2 , iy/ Xw> , where 
x V v' = (H 2 -!^') an d u vv i = lo v — io v i are the surface-atom 
dipole matrix element and the translational transition 
frequency, respectively. Hence, the coupling coefficient 
g vv i depends on the dipole matrix element x uv > and the 
transition frequency oj vu i. Since oj vv = 0, we have g vv = 
0. 

We note that the Hamiltonian of the ^-polarized acous- 
tic phonons is given by 



(10) 



The total Hamiltonian of the atom-phonon system is 

H = H A + H I + H B . (11) 

We use the above Hamiltonian to study the phonon- 
mcdiatcd decay of the atom. 

III. DYNAMICS OF THE ATOM 

In this section, we present the basic equations for the 
phonon-mediated decay processes. We derive a general 
master equation for the reduced density operator of the 
atom in subsection IIII A[ obtain analytical expressions 
for the relaxation rates and frequency shifts in subsec- 
tion IIII B[ and calculate the rates and the shifts in the 
framework of the Debye model in subsection IIII CI 

A. Master equation 

In the Heisenberg picture, the equation for the phonon 
operator i> q (i) is 

6 q (t) = -z^ q 6 q (i) - -^=5(t), (12) 

which has a solution of the form 

b q (t) = 6 q (t )e-^ ( *- to) - »W q (t). (13) 
Here, to is the initial time and Wq is given by 



Wq(i) 



-^-r)s( T )dT. (14) 



v / 2MNh ' 



Consider an arbitrary atomic operator O which acts only 
on the atomic states but not on the phonon states. The 
time evolution of this operator is governed by the Heisen- 
berg equation 

! ^ = L[H A {t)+H I (t),0{t)l (15) 
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which, with account of Eqs. ([6]) and (|13[) . yields 
dO(t) 



at =nl H *®>°™ 



+E 

q 

-E 



z [S(t),O(t)][b^(t )e 



-iuiqit—to) 



iW q (t)] 



:[6t(t ) e *--(*-*o)+iw3(t)][0(t),5(t)]. 



(16) 



We assume the initial density of the atom-phonon sys- 
tem to be the direct product state 



Pe(*o) = p{h)pB{h), 



(17) 



with the atom in an arbitrary state p(to) and the phonons 
in a thermal state 



p B {t ) = ^- 1 exphi/ B (t )/fc B T]. 



(18) 



Here, Z is the normalization constant and T is the tem- 
perature of the phonon bath. For the initial condition 
(|17p. the Bogolubov's lemma [Hj], applied to an arbitrary 
operator <d(t), asserts the following: 



<e(i)6 q (i )> = n q <[& q (to),e(t)]), 



(19) 



where the mean number of phonons in the mode q is 
given by 



1 



exp(hujq/ ksT) — 1 



(20) 



Let be an atomic operator. We then have the commu- 
tation relation [& q (i), ©(£)] = 0, which yields 



[6 q (i ) ) e(«)]=ie^(*- t °)[w q (i) ) e(i)]. 

Combining Eq. (19]) with Eq. (21]) leads to 



(e(t)& q (*o)> = ie^-^-'^nqdWqW.eW] 



(21) 



(22) 



We perform the quantum mechanical averaging for ex- 
pression (Til)]) and use Eq. (2"2"j) to eliminate the phonon 
operators & q (to) an d & q (io)- The resulting equation can 
be written as 



1 



[5(t),0(t)]W q (t)+W3(t)[0(t),5(t)]) 



E -^{w^mim + [s(t),o( t )]wi(t)). 



(23) 

We note that Eq. (23]) is exact. It does not contain 
phonon operators explicitly. The dependence on the 
phonon operators is hidden in the time shift of the oper- 
ator S(t) in expression (TJJ) for the operator Wq(t). 



We now show how the dependence of the operator 
Wq(t) on the phonon operators can be approximately 
eliminated. We assume that the atom-phonon coupling 
coefficients g vv i are small. The use of the zeroth-order 
approximation <j vv i(t) = <J 1/U i{t)e w " v '^ r ^ t ' in the expres- 
sion for S(t) [see Eq. 0] yields 



S(r) 



vv' 



(24) 



which is accurate to first order in the coupling coeffi- 
cients. Inserting Eq. ([24]) into Eq. (T4]) gives 



Wq(t) 



2tt 



^ 9vv> Oyv' (t)5- {ljJ v i v - UJg) , (25) 



where 



lim — 



e -i(u+ie)T dj . 



(26) 



i P 
2n u) 

Here, in order to take into account the effect of adiabatic 
turn-on of interaction, we have added a small positive 
parameter e to the integral and have used the limit to — > 
— oo. Introducing the notation 



W a 



2tt 



E 



guv'(Tuv'S-(Vu'u — w q ), (27) 



we can rewrite Eq. (23]) in the form 



(n q + l)<[S(t), 0(t)]X q (t) + KUt)p(t), S(t)}} 



+E^ 

q 

+ E ^(Kq{t)[0{t), S{t)] + [S(t), 0(t)]Kl(t)). (28) 
q 

In order to examine the time evolution of the reduced 
density operator pit) of the atom in the Schrodinger 
picture, we use the relation (0(t)) = Tr[O(t)p(0)] = 
Tr[O(0)p(t)], transform to arrange the operator O(0) at 
the first position in each operator product, and eliminate 
O(0). Then, we obtain the Liouville master equation 



dp® 
dt 



-^\H A ,p{t)] 



+ EK + mKqP(t),S] + [S, p(t)K\}} 
q 

+ J2 *«{[S, p(«)*q] + S]}- (29) 

q 

Equations ([28]) and ([29]) are valid to second order in 
the coupling coefficients. These equations allow us to 
study the time evolution and dynamical characteristics of 
the atom interacting with the thermal phonon bath. We 
note that Eq. (2"9"j) is a particular form of the Zwanzig's 
generalized master equation, which can be obtained by 
the projection operator method [l~3l |. 
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B. Relaxation rates and frequency shifts 



We use Eq. ([29)1 to derive an equation for the matrix 
elements pj.y = (j\p\j') of the reduced density operator 
of the atom. The result is 



dp. 



-ico 



- YM» + TiJP*' + + ^'M' ( 3 °) 



where the coefficients 



q 

5 + (oj, 



9jv9j'u'{S-(^ v j - W q ) 



•9jn9viiS- (u> utM - w q ) (31) 



and 



Jjj'uu' = 27r^- St 9jV3j' I /'[^-(^' W ' -w q ) 



+ - w q )] 

w, 



7£> = 27r X! ^r^M^M^^^ ~ W q) ( 32 ) 

are the decay parameters associated with the phonon 
emission and absorption, respectively. Here, the nota- 
tion d+(u) = 6*L(uj) has been used. 

Equation (|30[) describes phonon-induced variations in 
the populations and coherences of the translational levels 
of the atom. We analyze the characteristics of the relax- 
ation processes. For simplicity of mathematical treat- 
ment, we first consider only transitions from discrete lev- 
els. The equation for the diagonal matrix element p 3 j for 
a discrete level j can be written in the form 



dp 3 



Of 



(fii + ^ + c - c -)^- 

off-diagonal terms. 



(33) 



When the off-diagonal terms arc neglected, Eq. (|33|) re- 
duces to a simple rate equation. It is clear from Eq. (|33j) 
that the rate for the downward transition from an upper 
level I to a lower level k (k < I) is 



R ki = IkkU = 2lr 2^2 



nq + 1 



(34) 



while the rate for the upward transition from a lower level 
k to an upper level I (I > k) is 



Rtk = iukk = ^J2 rrMuik - (35) 



Equations (|34[) and (|35l) are in agreement with the re- 
sults of Gortel et at [10|], obtained by using the Fermi 
golden rule. We note that R kl and Rf k with I < k are 
mathematically equal to zero because they have no phys- 
ical meaning. For convenience, we introduce the notation 
Rik = Rf k , Rf k , or for I < k, I > k, or I = k, respec- 
tively. It is clear that the off-diagonal coefficients Rik 
with I 7^ k are the rates of transitions. However, the di- 
agonal coefficients Rkk have no physical meaning and are 
mathematically equal to zero. 

As seen from Eq. (|3"3"|) , the phonon- mediated depletion 
rate of a level k is Tkk — 2Re(7^ fc + l kk )- The explicit 
expression for this rate is 



l ~~9k^{^kii - Wq) 



1M 



q 2 xi \ 



(36) 



We note that T kk = E^ fe + R%k) 
write T kk = T e kk 1 r " 



■T% k) where 



kk 



H<k 



lik 



(37) 



and 



T^a \ ^ pa 

1 kk — V k 

M>fc 



(38) 



are the contributions due to downward transitions 
(phonon emission) and upward transitions (phonon ab- 
sorption), respectively. In the above equations, the sum- 
mation over p can be extended to cover not only the 
discrete levels but also the continuum levels. 

Meanwhile, the equation for the off-diagonal matrix 
element pik for a pair of discrete levels I and k can be 
written in the form dpik/dt = -{ioJik + "fa + 7$ + lt* k + 
~ftl)pik + . . . , or, equivalently, 



dp 



ik 



dt 



-i(tJik + A ik - iTik)pik + 



(39) 



Here the frequency shift Ai k is given by 



71q + 1 



9u 



tOl^j, — Wq LU^k + W q 



+ 



-I 2 
ipk 



Wq V w lv + w q w nk ~ W c 



(40) 



while the coherence decay rate Ti k is expressed as 

Tik = 7T 2~2, " q + 1 ^ w q) + 9kn s (vk» - w q )] 

q^ Wq 

+ ^E ~ [^(^ ~ ^ + 9lk 6 {^Hk - Wq)] ■ 

(41) 
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When we set I = k in Eq. (|40|). we find A kk = 0. 
When we set I = k in Eq. I|4ip. we recover Eq. (f3"6")) . 
We note that T ;fc - E M (^i + R% k + R a ,i + Rfr.)/2 = 
J^^iRfii +-R/i/fe)/2- Comparison between Eqs. and 
(l3"o) yields the relation Tife = (r« + Tkk)/2. We can also 
write T lk = Tf k + Tf k , where rf fc = + ^ fc )/2 

and r" fe = $Z„(i2£j + R^ k )/2 are the contributions due 
to downward transitions (phonon emission) and upward 
transitions (phonon absorption), respectively. In the 
above equations, the summation over p can be extended 
to cover not only the discrete levels but also the contin- 
uum levels. 

We now discuss phonon-mediated transitions from con- 
tinuum (free) levels. We start by considering free-to- 
bound transitions. For a continuum level / with energy 
£f > 0, the center-of-mass wave function <ff(x) is nor- 
malized per unit energy. In this case, the quantity R v f 
becomes the density of the transition rate. A free level / 
can be approximated by a level of a quasicontinuum 14 1 . 
A discretization of the continuum can be realized by using 
a large box of length L with reflecting boundary condi 



tions [15j. We label E n the energies of the eigenstates 
in the box and 4> n (x) the corresponding wave functions. 
Note that such states are standing-wave states 
The relation between a quasicontinuum-state wave func- 
tion <J) nf (x) , normalized to unity in the box, and the cor- 
responding continuum-state wave function (pf(x), nor- 
malized per unit energy, with equal energies E rif = £f, 
is Ql 



tp f {x) = 



dE n 



dn f 
L 



- -1/2 

4>n f (x) 
1/2 / \ 1/4 

j m 

\2E nf 



<t>n r {x). (42) 



Consequently, for a single atom initially prepared in the 
quasicontinuum standing- wave state |n/) = \4> nf ), the 
rate for the transition to an arbitrary bound state \v) is 
approximately given by 



Cr„/ - —VfH v f, 



(43) 



where vj = (2£f/m) 1 / 2 is the velocity of the atom in the 
initial standing- wave state |/). The phonon-mediated 
free-to-bound decay rate (adsorption rate) is then given 
by 



(44) 



where the summation includes only bound levels. It is 
clear from Eq. (|43[) that, in the continuum limit L — > oo, 
the rate G„/ tends to zero. This is because a free atom 
can be anywhere in free space and therefore the effect of 
phonons on a single free atom is negligible. 

In order to get deeper insight into the free-to-bound 
transition rate density R v f, we consider a macroscopic 



atomic ensemble in the thermodynamic limit [l4(. Sup- 
pose that there are Nq atoms in a volume with a large 
length L and a transverse cross section area So- Assume 
that all the atoms are in the same quasicontinuum state 
\nf) and interact with the dielectric independently The 
rate for the transitions of the atoms from the quasicon- 
tinuum state |n/) to an arbitrary bound state \v), defined 
as the time derivative of the number of atoms in the state 
\v), is D v f = N G v f. In order to get the rate for the con- 
tinuum state |/), we need to take the thermodynamical 
limit, where L — * oo and Nq — > oo but Nq/L remains 
constant. Then, the rate for the transitions of the atoms 
from the continuum state |/) to an arbitrary bound 
state \v) is given by D v j — nhpoSoVfR^f — 2irhAl fRuf ■ 
Here, pa = Nq/LSq is the atomic number density and 
Mf = poSoVf /2 is the number of atoms incident into the 
dielectric surface per unit time. It is clear that the tran- 
sition rate is proportional to the incidence rate A// 
as well as the transition rate density R v f. We emphasize 
that -D„/ is a characteristics for a macroscopic atomic en- 
semble in the thermodynamic limit while G„/ is a mea- 
sure for a single atom. When the length of the box, L, 
and the number of atoms, Nq, are finite, the dynamics of 
the atoms cannot be described by the free-to-bound rate 
D v f directly. Instead, we must use the transition rate 
per atom G v j — D v f/No, which depends on the length 
L of the box that contains the free atoms [see Eq. (|4"3"|) ]. 

In a thermal gas, the atoms have different velocities 
and, therefore, different energies. For a thermal Maxwell- 
Boltzmann gas with temperature To, the distribution of 
the kinetic energy £ f of the atomic center-of-mass motion 
along the x direction is 



P(£f) 



i 



-£f/k B T 



Vr7kB% \f£ 



J 



(45) 



The transition rate to an arbitrary bound state \v) is 
then given by G v t = J °° G„fP{£f) d£f, i.e. 



G vTo = ^-fe- £ t^ T °R uJ d£ f 



(46) 



where Ad = {2ttTl 2 /mkBTo) 1 / 2 is the thermal de Broglie 
wavelength. The phonon-mediated free-to-bound decay 
rate (adsorption rate) is given by 



/>oo 

G To = ^G, To = / G f P(£ f ) 
Jo 



d£ 



(47) 



In the above equation, the summation over v includes 
only bound levels. Note that Eq. (|46l is in qualitative 
agreement with the results of Refs. [a. X\M±. 

It is easy to extend the above results to the case of 
free-to- free transitions. Indeed, it can be shown that the 
density of the rate for the transition from a quasicontin- 
uum state \rif), which corresponds to a free state |/), to 
a different free state |/') is given by 



Qff 



7i h 



- v f R f'f- 



(48) 
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For convenience, we introduce the notation Qff = Qff 
or for £f> < £f or £f> > £f, respectively, and QJ/f — 
Qff or for £fi > £f oy £f < £/, respectively. Then, we 
have Qff — Qff, 0, or Q°f, f for £f < £f, £f — £f, or 
£f > £f, respectively. The downward (phonon-emission) 
and upward (phonon-absorption) free-to-free decay rates 
for the free state |/) are given by 



Q), s d£f 



(49) 



over a sphere of radius qo = {Qit 2 N /V) 1 / 3, , where V is the 
volume of the solid. The Debye frequency and the Debye 
temperature are given by ujjj = vqu and Tjj = TiuJo/kB, 
respectively. For fused silica, we have v = 5.96 km/s, 
NM/V = 2.2 g/cm 3 , and M = 9.98 x 1(T 26 kg 0. Us- 
ing these parameters, we find qn — 109.29 x 10 6 cm -1 , 
uj d = 10.4 THz, and T D = 498 K. In order to perform 
the summation over phonon states in the framework of 
the Debye model, we invoke the thermodynamic limit, 
i.e., replace 



and 



Qffd£f, 



(50) 



respectively. The total free-to-free decay rate for the free 
state I/) is Q f = Q e f + Q f = f™Q f , f d£ r . 

For a thermal gas, we need to replace the transition 
rate density Qff and the decay rate Qf by QfT = 
Jp°° Qf'f p ( £ f) d£ f and Qt = /(T 'QfP(£f)d£ f , respec- 
tively, which are the averages of Qff and Q / , respec- 
tively, with respect to the energy distribution P(£f) 
of the initial state. Like in the other cases, we have 

Q/'To = Q/'To + QfT and QTo = Qt + Qt > where 

Q/'To = / QffP(£f)d£f, 



if To 



Q a flf P{£f)d£ f 



(51) 



are the downward and upward transition rate densities 
and 



Qt 



Qt„ 



Q}P(£ f )d£ f , 
Q a f P(£f)d£f 



(52) 



are the downward and upward decay rates. The thermal 
decay rates Qt an d Qt describe the cooling and heating 
processes, respectively. It can be easily shown that Q^ < 
Qt , Qt > Qt , and Q e To = Q a TQ when T < T, T ° > 
T, and T = T, respectively. The relation Q T < Q^ 
(Qt > Qt )> obtained for T < T (T > T), indicates 
the dominance of heating (cooling) of free atoms by the 
surface. 



C. Relaxation rates and frequency shifts in the 
framework of the Debye model 

In order to get insight into the relaxation rates and 
frequency shifts, we approximate them using the Debye 
model for phonons. In this model, the phonon frequency 
(jj q is related to the phonon wave number q as o> q = v q, 
where v is the sound velocity. Furthermore, the summa- 
tion over the hrst Brillouin zone is replaced by an integral 



V 

8^ 



, dq — 



q|<9-o 



3iV 



. Uj'idldc, 



(53) 



Then, for transitions between an upper level I and a lower 
level fc, where < % < lj d , Eqs. (fM)) and ([33]) yield 



R 



37T 9 



and 



no 



3tt 



lk MTiu> 3 D 



ci2 



(54) 



(55) 



Here, nik is given by Eq. ([2H)l with w q replaced by u>ik. 
We emphasize that, according to Eqs. (|54|) and (|55|l . 
the phonon-emission rate and the phonon-absorption 
rate Rf k depend not only on the matrix element -F^ of 
the force but also on the translational transition fre- 
quency ujik- The frequency dependences of the transi- 
tion rates are comprised of the frequency dependences 
of the mean phonon number n/^ , the phonon mode den- 
sity 3Nu>f k /oji), and the matrix element Fi 



lk 



—mujf k xik of the force. An additional factor comes from 
the presence of the phonon frequency in Eq. ([5]) for 
the surface displacement and, consequently, in the atom- 
phonon interaction Hamiltonian ([5]). It is clear that an 
increase in the phonon frequency leads to a decrease in 
the mean phonon number and an increase in the phonon 
mode density. The matrix element of the force usu- 
ally first increases and then decreases with increasing 
phonon frequency. Due to the existence of several com- 
peting factors, the frequency dependences of the tran- 
sition rates are rather complicated. They usually first 
increase and then decrease with increasing phonon fre- 
quency. We note that, for transitions with u>ik > oju, we 



have Ri 



Da 
H lk 



0. 



We conclude this section by noting that the use of Eq. 
53"| in Eq. (|4TI)) yields the frequency shift 



AO) 

*lk 



AT) 
*lk > 



where 



A (o) _ 



E 



F 2 

Ulu - LO 



F 2 

1 ^k 



(56) 



jduj (57) 
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and 

A (T) 3 \ - f ( u^F t u^kF k \_ 

M 

(58) 

are the zero- and finite-temperature contributions, re- 
spectively. In Eq. (|58|) . n w is given by Eq. ([20]) with w q 
replaced by u). 



IV. NUMERICAL RESULTS AND 
DISCUSSIONS 

In this section, we present the numerical results based 
on the analytical expressions derived in the previous 
section for the phonon-mediated relaxation rates of the 
translational levels of the atom. In particular, we use 
Eqs. (fS"4")) and (|55[) . obtained in the framework of the De- 
bye model, for our numerical calculations. We consider 
transitions from bound states as well as free states. The 
transitions from bound states to other translational lev- 
els occur in the case where the atom is initially already 
adsorbed or trapped near the surface. The transitions 
from free states to other translational levels occur in the 
processes of adsorbing, heating, and cooling of free atoms 
by the surface. Due to the difference in physics of the ini- 
tial situations, we study the transitions from bound and 
free states separately. 

A. Transitions from bound states 



FIG. 3: Phonon-emission rates R1i v from the vibrational lev- 
els (a) v — 280 and (b) v = 120 to other levels v' as functions 
of the lower-level energy £„/ . The arrows mark the initial 
states. The parameters of the solid are M = 9.98 x 10~ 26 kg 
and u>d = 10.4 THz. The temperature of the phonon bath is 
T = 300 K. Other parameters are as in Fig. [1] 



FIG. 4: Phonon-absorption rates from the vibrational 

levels (a) v = 280 and (b) v = 120 to other levels v' as func- 
tions of the upper-level energy £„/ . The left (right) panel 
in each row corresponds to bound-to-bound (bound-to-free) 
transitions. The arrows mark the initial states. The param- 
eters used are as in Fig. [3] The temperature of the phonon 
bath is T = 300 K. 

We start from a given bound level and calculate the 
rates of phonon-mediated atomic transitions, both down- 
ward and upward. The profiles of the phonon-emission 
(downward-transition) rate R%i v [see Eq. (|54"|) ] and the 
phonon-absorption (upward-transition) rate R%i v [see Eq. 
(I55[) ] are shown in Figs. [3] and 21 respectively. The upper 



(lower) part of each of these figures corresponds to the 
case of the initial level v — 280 (v = 120), with energy 
£ v = -156 MHz [E v = -8.4 THz). The left (right) panel 
of Fig. 0] corresponds to bound-to-bound (bound-to-free) 
upward transitions. The temperature of the surface is 
assumed to be T = 300 K. As seen from Figs. [3] and 2J 
the transition rates have pronounced localized profiles. 
Due to the competing effects of the mean phonon num- 
ber, the phonon mode density, and the matrix element 
of the force, the transition rates usually first increase 
and then decrease with increasing phonon frequency. It 
is clear from a comparison of Figs. [3fa) and [3{b) and 
also a comparison of Figs.[4ja) and[4jb) that transitions 
from shallow levels have probabilities orders of magni- 
tude lower than those from deeper levels. The main rea- 
son is that the wave functions of the shallow states are 
spread further away from the surface than those for the 
deep states. Due to this difference, the effects of the sur- 
face vibrations are weaker for the shallow levels than for 
the deep levels. Another pertinent feature that should 
be noted from the figure is the following: Since transi- 
tion frequencies involved are large, they may overshoot 
the Debye frequency lod = 10.4 THz, leading to a cutoff 
on the lower (higher) side of the frequency axis for the 
emission (absorption) curve. 

In order to see the overall effect of the individual tran- 
sition rates shown above, we add them up. First we ex- 
amine the phonon-absorption rates of bound levels. The 
total phonon-absorption rate of a bound level v is the 
sum of the individual absorption rates R® over all the 
upper levels fi, both bound and free [see Eq. l[3"8"])]. We 
plot in Fig. [5] the contributions to T" v from two types 
of transitions, bound-to-bound and bound-to-free (des- 
orption) transitions. The solid curve of the figure shows 
that the bound-to-bound phonon-absorption rate is large 
(above 10 10 s _1 ) for deep and intermediate levels. How- 
ever, it reduces dramatically with increasing v in the 
region of large v and becomes very small (below 10~ 5 
s _1 ) for shallow levels. Meanwhile, the dashed curve of 
Fig. [5] shows that the bound-to-free phonon-absorption 
rate (i.e., the desorption rate) is zero for deep levels, since 
the energy required for the transition is greater than the 
Debye energy However, the desorption rate is sub- 
stantial (above 10 5 s _1 ) for intermediate and shallow lev- 
els. Thus, the total phonon-absorption rate T® u is mainly 
determined by the bound-to-bound transitions in the case 
of deep levels and by the bound-to-free transitions in the 
case of shallow levels. One of the reasons for the dramatic 
reduction of the bound-to-bound phonon-absorption rate 
in the region of shallow levels is that the number of up- 
per bound levels fx becomes small. The second reason is 
that the frequency of each individual transition becomes 
small, leading to a decrease of the phonon mode den- 
sity. The third reason is that the center-of-mass wave 
functions of shallow levels are spread far away from the 
surface, leading to a reduction of the effect of phonons 
on the atom. 

Unlike the bound-to-bound phonon-absorption rate, 
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the bound-to-free phonon-absorption rate is substantial 
in the region of shallow levels. This is because the free- 
state spectrum is continuous and the range of the bound- 
to-free transition frequency can be large (up to the De- 
bye frequency ujd — 10.4 THz). The gradual reduction of 
the bound-to- free phonon-absorption rate in the region of 
shallow levels is mainly due to the reduction of the time 
that the atom spends in the proximity of the surface. 



transitions, which correspond to the adsorption process. 
According to Eq. (j4*3"j) . the free-to-bound (more exactly, 
quasicontinuum-to-bound) transition rate G„/ depends 
not only on the continuum-to-bound transition rate den- 
sity R v f but also on the length L of the free-atom quan- 
tization box. To be specific, we use in our numerical 
calculations the value L = 1 mm, which is a typical size 
of atomic clouds in magneto-optical traps [17[ • 



FIG. 5: Contributions of bound-to-bound (solid curve) and 
bound-to-free (dashed curve) transitions to the total phonon- 
absorption rate r°„ versus the vibrational quantum number 
v of the initial level. The parameters used are as in Fig. [3] 
The temperature of the phonon bath is T = 300 K. 

The total phonon-emission rate Y e uv [see Eq. ([3"T)l ] and 
the total phonon-absorption rate [see Eq. (|38|) ] are 
shown in Fig. [5] by the solid and dashed curves, respec- 
tively, ft is clear from the figure that emission is com- 
parable to but slightly stronger than absorption. Such a 
dominance is due to the fact that phonon emission moves 
the atom to a center-of-mass state closer to the surface 
while phonon absorption changes the atomic state in the 
opposite direction (see Figs. [1] and [2]). Our results for the 
rates are in good qualitative agreement with the results 
of Oria et ai, albeit with the Morse potential [5(. We 
stress that we include a large number of vibrational lev- 
els as a consequence of the deep silica-cesium potential. 
Note that the earlier work on this theme involved much 
fewer levels [5[. 



FIG. 6: Phonon-emission decay rate FJ„ (solid lines) and 
phonon-absorption decay rate r"„ (dashed lines) of a bound 
level as functions of the vibrational quantum number v. The 
inset shows the rates in the linear scale to highlight the dif- 
ferences in the dissociation limit. The parameters used are as 
in Fig. [3] The temperature of the phonon bath is T = 300 K. 



FIG. 7: Same as in Fig. |6] except that T = 30 K. 

We next study the effect of temperature on the decay 
rates. The results for the phonon-mediated decay rates 
for T = 30 K are shown in Fig. [7] In contrast to Fig. [51 
the absorption rate is now much smaller than the corre- 
sponding emission rate for both shallow and deep levels. 
Thus, while it is difficult to distinguish the two log-scale 
curves for deep and shallow levels at room temperature 
(see Fig. [5J), they are well resolved at low temperature. 

B. Transitions from free states 

We now calculate the rates for transitions from free 
states to other levels. We first examine free-to-bound 



FIG. 8: Free-to-bound transition rates G v f for transitions 
from the free plane-wave states with energies (a) Bf = 2 MHz 
and (b) Bf = 3.1 THz to bound levels v as functions of the 
bound-level energy £„. The arrows mark the energies of the 
initial free states. The insets show G„/ on the log scale versus 
£ v in the range from -200 MHz to -0.2 MHz to highlight 
the rates to shallow bound levels. The length of the free- 
atom quantization box is L = 1 mm. The temperature of the 
phonon bath is T — 300 K. Other parameters are as in Fig. [3] 

We plot in Fig. [5] the free-to-bound transition rate G„/ 
[see Eq. (|43[) ] as a function of the vibrational quantum 
number v. The upper (lower) part of the figure corre- 
sponds to the case of the initial-state energy Ef = 2 
MHz (Ef = 3.1 THz), which is close to the average ki- 
netic energy per atom in an ideal gas with temperature 
T = 200 iaK (T = 300 K). We observe that the free-to- 
bound transition rate first increases and then decreases 
with increasing transition frequency ujfy = (Ef — E v )/h. 
Such behavior results from the competing effects of the 
mean phonon number, the phonon mode density, and the 
matrix element of the force, like in the case of bound- 
to-bound transitions (see Fig. [3]). We also see a cut- 
off of the transition frequency, which is associated with 
the Debye frequency. Comparison of Figs. Ufa) and[5fb) 
shows that the transitions from low-energy free states 
have probabilities orders of magnitude smaller than those 
from high-energy free states. One of the reasons is that 
the transition rate G v f is proportional to the velocity 
v f = (2Ef/mf/ 2 [sec Eq. (@3J)]. The dependence of 
the transition rate density i?„/ on the transition fre- 
quency ujf u also plays an important role. Because of this, 
the rates for the transitions from low-energy free states 
to shallow bound levels are very small [see the inset of 
Fig. Ha)]. 



FIG. 9: Frec-to-bound decay rate G/ as a function of the 
free-state energy Bf. The inset highlights the magnitude and 
profile of the decay rate for Bf in the range from to 20 
MHz. The temperature of the phonon bath is T — 300 K. 
Other parameters are as in Fig. [5] 

We show in Fig. [9] the free-to-bound decay rate G/ 
[see Eq. (|44|>] . which is a characteristic of the adsorp- 
tion process, as a function of the free-state energy Ef. 
We see that G/ first increases and then decreases with 
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increasing £f. The increase of Gf with increasing £f in 
the region of small £ t (see the inset) is mainly due to 
the increase in the atomic incidence velocity Vf. In this 
region, we have Gf oc Vf oc ^J~£~f [see Eqs. (|4l?)) and 
flS}]. For £/ in the range from to 20 MHz, which is 
typical for atoms in magneto-optical traps, the maximum 
value of Gf is on the order of 10 4 s _1 (see the inset of 
Fig. [5]). Such free-to-bound (adsorption) rates are sev- 
eral orders of magnitude smaller than the bound-to-free 
(desorption) rates (see the dashed curve in Fig. [5]). The 
decrease of Gf with increasing £f in the region of large 
£f is mainly due to the reduction of the atom-phonon 
coupling coefficients. 



FIG. 10: Free-to-bound transition rates G„t for transitions 
from the thermal states with temperatures (a) To = 200 /iK 
and (b) To = 300 K to bound levels v as functions of the 
bound-level energy £ „ . The insets show G„t on the log scale 
versus £ v in the range from —200 MHz to —0.2 MHz to high- 
light the rates to shallow bound levels. The temperature of 
the phonon bath is T — 300 K. Other parameters are as in 
Fig.il 



FIG. 11: Free-to-bound decay rate Gt as a function of the 
atomic temperature To in the ranges (a) from 100 /iK to 400 
(iK and (b) from 50 K to 350 K. The temperature of the 
phonon bath is T = 300 K. Other parameters are as in Fig. [8] 

In a thermal gas, the adsorption process is charac- 
terized by the transition rate G„t [see Eq. (f4"o]) ] and 
the decay rate Gt [see Eq. (jT7|) ]. which are the av- 
erages of the free-to-bound transition rate G„/ and the 
free-to-bound decay rate G/, respectively, over the free- 
state energy distribution (|4"5")) . We plot the free-to-bound 
transition rate G v t q and the free-to-bound decay rate 
Gt in Figs. [TU] and QTJ respectively. Comparison be- 
tween Figs-flUTa - ) and^a) shows that the transition rates 
from low-temperature thermal states and low-energy free 
states look quite similar to each other. The reason is that 
the spread of the energy distribution is not substantial 
in the case of low temperatures. The spread of the en- 
ergy distribution is however substantial in the case of 
high temperatures, leading to the softening of the cut- 
off frequency effect [compare Fig. HOTb') with Fig. [9fb)]. 
Figure QT] shows that the free-to-bound decay rate Gt 
first increases and then reduces with increasing atomic 
temperature To. For To in the range from 100 fjK to 400 
/iK, which is typical for atoms in magneto-optical traps, 
the maximum value of Gt is on the order of 10 4 s _1 [see 
Fig. fTlTa)]. Such free-to-bound (adsorption) rates are 
several orders of magnitude smaller than the bound-to- 
free (desorption) rates (see the dashed curve in Fig. [5]). 
Figure [TlT a) shows that, in the region of low atomic tem- 
perature T , one has Gt oc \fT$, in agreement with the 
asymptotic behavior of Eqs. (|4"o| and (ITT)) . 



FIG. 12: Free-to- free transition rate densities Qff for the 
upward (solid lines) and downward (dashed lines) transitions 
from the free states |/) with energies (a) £f = 2 MHz and (b) 
£f — 3.1 THz to other free states |/') as functions of the final- 
level energy £ii. The arrows mark the energies of the initial 
free states. The inset in part (a) shows Qff versus Eti in 
the range from to 4 MHz to highlight the small magnitude 
of the rate density for downward transitions (dashed line). 
The temperature of the phonon bath is T = 300 K. Other 
parameters are as in Fig. [5] 

We now examine free-to-free transitions, both upward 
and downward, which corresponding to the heating and 
cooling processes of free atoms by the surface. We plot in 
Fig. HH the free-to-free transition rate density Qff [see 
Eq. (|48)) ] as a function of the final-level energy £f>. The 
upper (lower) part of the figure corresponds to the case 
of the initial-state energy £f — 2 MHz (£f = 3.1 THz), 
which is close to the average kinetic energy per atom in 
an ideal gas with temperature T = 200 [iK (T = 300 
K). The rate densities are shown for the upward (phonon- 
absorption) and downward (phonon-emission) transitions 
by the solid and dashed lines, respectively. The fig- 
ure shows that the free-to-free transition rate density in- 
creases or decreases with increasing transition frequency 
if the latter is not too large or is large enough, respec- 
tively. We also observe a signature of the Debye cutoff 
of the phonon frequency. Comparison of Figs. I12f a) and 
|T2f b) shows that transitions from low-energy free states 
have probabilities orders of magnitude smaller than those 
from high-energy free states. Figure [12^ a) and its inset 
show that, when the energy of the free state is low, the 
free-to-free downward (cooling) transition rate is very 
small as compared to the free-to-free upward (heating) 
transition rate. 



FIG. 13: Free-to-free upward and downward decay rates Q°j 
(solid lines) and Qf (dashed lines) as functions of the energy 
£f of the initial free state. The insets highlight the magni- 
tudes and profiles of the decay rates for £f in the range from 
to 20 MHz. The temperature of the phonon bath is T = 300 
K. Other parameters are as in Fig. [8] 

We show in Fig. [T2] the free-to-free upward (phonon- 
absorption) and downward (phonon-emission) decay 
rates Q a f [see Eq. ((50)) ] and Q e j [see Eq. ([4Tj]) ] as functions 
of the free-state energy £f. We observe that QJ and Q/ 
increase with increasing £f in the range from to 8 THz. 
The increase of Q a j with increasing £f in the region of 
small £ f (see the left inset) is mainly due to the increase 
in the atomic incidence velocity Vf. In this region, we 
have Q a f oc v f oc y/£j [see Eqs. gH]) and ([50])]. The 
increase of Qff with increasing £f in the region of small 
£f (see the right inset) is due to not only the increase in 
the atomic incidence velocity Vf [see Eq. (|4"5)) ] but also 
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the increase of the transition rate density Q e fif and the 
increase of the integration interval (0, £f) [see Eq. (|4"9j) ]. 
In this region, the dependence of on the energy £f 

is of higher order than ■ The left inset of Fig. [T3l 
shows that, for £f in the range from to 20 MHz, the 
maximum value of QJ is on the order of 10 4 s _1 . Such 
free-to-free upward (heating) decay rates are comparable 
to but about two times smaller than the corresponding 
free-to-bound (adsorption) decay rates (see the inset of 
Fig. [9]). Meanwhile, the right inset of Fig. [T3l shows that, 
in the region of small £/, the free-to- free downward (cool- 
ing) decay rate Q % is very small. 



FIG. 14: Free-to-free transition rate densities Q}t for up- 
ward transitions (solid lines) and Q}t f° r downward tran- 
sitions (dashed lines) from the thermal states with tempera- 
tures (a) T = 200 /jK and (b) T = 300 K to free levels / 
as functions of the free-level energy £f. The inset in part (a) 
shows the rate densities versus £f in the range from to 8 
MHz to highlight the small magnitude of Q/t (dashed line). 
The temperature of the phonon bath is T = 300 K. Other 
parameters are as in Fig. [8] 



line), while the cooling rate Qy o (dashed line) is negligi- 
ble. Figure [T5Ta) shows that, in the region of low atomic 
temperatures, one has Qt — Qt ^ V^o? m agreement 
with the asymptotic behavior of expressions ([52")) . The 
figure also shows that Qj, quickly increases with increas- 
ing atomic temperature To. The relation Qf Ta < Q^ , 
obtained for To < T, indicates the dominance of heating 
of cold free atoms by the surface. The substantial mag- 
nitude of the free-to-bound transition rate Gt (dotted 
line) indicates that a significant number of atoms can be 
adsorbed by the surface. According to Fig. [TBT b). the 
free-to-free downward transition rate Q^ o (dashed line) 
crosses the upward transition rate Q^ Q (solid line) when 
T = T = 300 K, and then becomes the dominant decay 
rate. The relation Q^ > Qt , obtained for T > T, indi- 
cates the dominance of cooling of hot free atoms by the 
surface. 



V. CONCLUSIONS 



FIG. 15: Free-to-free decay rates Qt (solid lines) and 
Qt (dashed lines) for upward and downward transitions, re- 
spectively, as functions of the atomic temperature To in the 
ranges (a) from 100 fjK to 400 fiK and (b) from 50 K to 350 
K. For comparison, the free-to-bound decay rate Gt is re- 
plotted from Fig. [11] by the dotted lines. The temperature of 
the phonon bath is T — 300 K. Other parameters are as in 
Fig.® 

In the case of a thermal gas, the phonon-mediated heat 
transfer between the gas and the surface is characterized 
by the free-to-free transition rate densities Qjt an d Q/t 
[see Eqs. ([5"Tj) ] and the free-to- free decay rates Qy o and 
Qt [ see -^qs. ([52]) ]. We plot the free-to-free transition 
rate densities QfT an d Qm m Fig- HU Comparison 
between Figs. [HT a) and[T27a) shows that the transition 
rate densities from low-temperature thermal states and 
low-energy free states are quite similar to each other. 
The spread of the initial-state energy distribution is not 
substantial in this case. However, the energy spread of 
the initial state is substantial in the case of high tem- 
peratures, concealing the cutoff frequency effect [com- 
pare Fig. [TJJb) with Fig. H5Tb)]. We display the free- 
to-free decay rates Q^ g and Qf^ in Fig. [15] The solid 
and dashed lines correspond to the upward (heating) and 
downward (cooling) transitions, respectively. For com- 
parison, the free-to-bound decay rate (adsorption rate) 
Gt is re-plotted from Fig. [TT] by the dotted lines. We 
observe that, for T in the range from 100 /iK to 400 /iK 
[see Fig. ITBT a)]. the adsorption rate Gt (dotted line) is 
about two times larger than the heating rate Q?p (solid 



In conclusion, we have studied the phonon-mediated 
transitions of an atom in a surface-induced potential. 
We developed a general formalism, which is applicable 
for any surface-atom potential. A systematic derivation 
of the corresponding density-matrix equation enables us 
to investigate the dynamics of both diagonal and off- 
diagonal elements. We included a large number of vi- 
brational levels originating from the deep silica-cesium 
potential. We calculated the transition and decay rates 
from both bound and free levels. We found that the 
rates of phonon-mediated transitions between transla- 
tional levels depend on the mean phonon number, the 
phonon mode density, and the matrix element of the force 
from the surface upon the atom. Due to the effects of 
these competing factors, the transition rates usually first 
increase and then reduce with increasing transition fre- 
quency. We focused on the transitions from bound states. 
Two specific examples, namely, when the initial level is 
a shallow level also when it can be one of the deep levels 
have been worked out. We have shown that there can be 
marked differences in the absorption and emission behav- 
ior in the two cases. For example, both the absorption 
and emission rates from the deep bound levels can be sev- 
eral orders (in our case, six orders) of magnitude larger 
than the corresponding rates from the shallow bound lev- 
els. We also analyzed various types of transitions from 
free states. We have shown that, for thermal atomic ce- 
sium with temperature in the range from 100 /iK to 400 
/iK in the vicinity of a silica surface with temperature of 
300 K, the adsorption (free-to-bound decay) rate is about 
two times larger than the heating (free-to- free upward de- 
cay) rate, while the cooling (free-to-free downward decay) 
rate is negligible. 
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